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Isoelectronic differential force measurements provide a unique opportunity to probe controversial 
features of the thermal Casimir effect, that are still much debated in the current literature. Isolec- 
tronic setups offer two major advantages over conventional Casimir setups. On one hand they are 
immune from electrostatic forces caused by potential patches on the plates surfaces, that plague 
present Casimir experiments especially for separations in the micron range. On the other hand they 
can strongly enhance the discrepancy between alternative theoretical models that have been pro¬ 
posed to estimate the thermal Casimir force for metallic and magnetic surfaces. Thanks to these two 
features, isoelectronic differential experiments should allow to establish conclusively which among 
these models correctly describes the thermal Casimir force. 
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I. INTRODUCTION 

Over sixty years ago [1] the dutch physicist Hendrik 
Casimir, building on the basic principles of Quantum 
Electrodynamics, predicted that two discharged perfectly 
conducting plane-parallel surfaces in vacuum attract each 
other with a (unit-area) force of magnitude 


This force orginates from zero-point quantum fluctua¬ 
tions of the electromagnetic (em) field that according 
to Quantum Theory fill empty space, or more precisely 
from modifications of the spectrum of these fluctuations 
caused by the presence of polarizable (but otherwise neu¬ 
tral) material surfaces. In his pioneering paper, Casimir 
considered two perfectly conducting plates at zero tem¬ 
perature. The theory of the Casimir effect for real mate¬ 
rial surfaces was developed a few years later by Lifshitz 
[2], by extending to dispersion forces Rytov’s theory of 
electromagnetic fluctuations [3]. For a general overview 
on the Casimir effect see W- 

During the fifty years following Casimir’s seminal pa¬ 
per, a few experiments were performed to observe the 
Casimir force which, apart from providing a qualitative 
confirmation of the effect, had the important merit of 
identifying the experimental problems that had to be ad¬ 
dressed for a successful observation of the tiny Casimir 
force. The modern era of the Casimir effect started in 
1997 with a landmark torsion-balance experiment by S. 
Lamoreaux [B], soon followed by the AFM experiment 
of Mohideen and Roy [S] ■ These experiments opened the 
era of precision Casimir experiments, by which it became 
possible for the first time to explore a number of features 
of the Casimir effect relating to material properties of 
the involved surfaces, that could be predicted on the ba¬ 
sis of Lifshitz theory. Several other experiments followed 
in rapid succession, which utilized both metallic surfaces 


[IOHM] , as well as surfaces made of diverse materials like 
semiconductors conductive oxides m, ferromag¬ 
netic metals HZi ng, and surfaces immersed in liquids 
[19] . Superconducting Casimir devices have been studied 
as well pniES] . Much excitement was spurred by the pos¬ 
sible exploitation of the Casimir force in the actuation of 
micro and nano machined devices [24) . which stimulated 
investigations of the Casimir effect for microstructured 
surfaces p51130j . For a review of these and many other 
experiments see miiTiin]- 

The specific subject of this paper is the effect on the 
Casimir force of the finite temperature T of the plates, 
a problem that has attracted a good deal of interest in 
recent years. In his seminal paper, Lifshitz showed that 
thermal fluctuations of the em field provide an extra con¬ 
tribution to the Casimir force, called thermal Casimir 
force. Surprisingly, estimating the magnitude of the ther¬ 
mal force for conducting plates turned into an intriguing 
puzzle. In essence the puzzle is about the role played 
by relaxation properties of conduction electrons in Lif¬ 
shitz theory. It turns out that Lifshitz formula predicts 
significantly different magnitudes for the thermal force 
depending on whether the optical data of the conductor 
are extrapolated towards zero frequency on the basis of 
the Drude model (which does take dissipation into ac¬ 
count) or instead by the dissipationless plasma model 
of IR optics. In addition to predicting different magni¬ 
tudes for the thermal force, it has been shown that the 
two prescriptions have important thermodynamic conse¬ 
quences: while the Drude prescription leads to a violation 
of Nernst heat theorem (in the idealized case of two con¬ 
ducting plates with a perfect crystal structure) [311133] j 
the plasma prescription violates the Bohr-van Leeuwen 
theorem of classical statistical physics [SU |35| . 

The experimental situation is contradictory. Several 
small distance experiments [TTHl4l 118] , probing separa¬ 
tions below one micron, appear to be in agreement with 
the plasma model, and to rule out the larger thermal 
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force predicted by the Drude model. These experiments 
provide the most precise measurements of the Casimir 
force to date, with errors in the percent range. One has 
to bear in mind however that for separations a smaller 
than the thermal length At = fic/{2'KkBT) (for T = 300 
K A 300 K = 1-2 /im) the thermal force only represents 
a small correction to the zero-point force, and therefore 
the theoretical interpretation of these experiments is very 
delicate. For a review of these experiments see |36] . 

In principle, observing the thermal Casimir force 
should be easier for separations a > At because for these 
large separations the thermal force is dominant. Discrim¬ 
inating between the Drude and the plasma models should 
be easier as well, because for large distances the two pre¬ 
scriptions predict markedly different magnitudes for the 
Casimir force. For example, for two plane-parallel con¬ 
ducting surfaces the Drude model predicts for a S> At a 
Casimir pressure of magnitude ({3)kBT/8Tra^, while the 
plasma model predicts a magnitude twice as large. Un¬ 
fortunately observing the thermal force for separations in 
the micron region is very difficult too, not only because 
the Casimir force quickly gets very small as the sepa¬ 
ration increases, but also because of unavoidable stray 
electrostatic forces that mask the Casimir force. These 
stray forces, that cannot be eliminated by applying a 
bias potential, originate from patches of varying poten¬ 
tial on the surfaces caused by spatial changes of crys¬ 
talline structure and/or by adsorbed impurities. Stray 
electrostatic forces over 100 times stronger than the ther¬ 
mal Casimir force were indeed observed in an experiment 
with A1 surfaces[37] in the range from 3.5 to 5 /rm. Large 
electrostatic forces were reported as well in a recent ex¬ 
periment by the Yale group [5B], which claims to have 
observed the thermal force between a large sphere and a 
plate both covered with gold, in the wide range of sep¬ 
arations form 0.7 to 7.3 /rm. The results have been in¬ 
terpreted by the authors as being in accordance with the 
Drude prescription. This experiment has been criticized 
[36j , because the thermal Casimir force was obtained only 
after subtracting from the total measured force the much 
larger electrostatic force. The subtraction was perfomed 
by making a fit of the total observed force, based on a 
two-parameter model of the electrostatic force, and not 
by a direct and independent measurement, as it would 
have been desirable. The problem of patch potentials is 
regarded as a major obstacle for present Casimir exper¬ 
iments, and dedicated techniques based on Kelvin probe 
force microscopy are being developed to achieve a direct 
observation of the patches with the necessary spatial res¬ 
olution [39]. 

The contradictory results of recent experiments call for 
new experiments specifically designed to probe the ther¬ 
mal Casimir force. Recently the author proposed two 
setups gnmi] that should allow for an unambiguous ob¬ 
servation of the thermal Casimir force. Both setups are 
based on isoelectronic differential force measurements, 
an idea pioneered by the lUPUI group in searches for 
non-newtonian gravitational forces in the sub-micron re¬ 


gion |3H US] . The unique virtue of this approach is that 
it is immune from the problem of patch potentials that 
plague conventional Casimir experiments, especially for 
large separations. A thorough analysis of the limitations 
on the sensitivity of isolelectronic experiments resulting 
from random spatial variations of patch forces has been 
recently carried out [44] . confirming the high suppres¬ 
sion of patch forces in isoelectronic setups. In iiniiii] 
it was proved that the isolelectronic technique provides 
a powerful tool to observe the elusive thermal Casimir 
force. In this paper we further develop the findings of 
gnmi], and we demonstrate that by suitably choosing 
the materials of the samples it is possible to strongly en¬ 
hance the discrepancy between the Drude and the plasma 
models, both for non-magnetic and for magnetic metals. 
This makes one confident that by this technique it should 
be possible to establish conclusively which among these 
models correctly describes the thermal Casimir force. In 
particular, it should be possible to clarify if and to what 
extent the large magnetic permeability of ferromagnetic 
materials influences the Casimir force, a problem that 
has been investigated recently by the Riverside group 
[mnH]. For ferromagnetic metals the isolectronic scheme 
is especially effective, because with such a setup the rel¬ 
ative difference between alternative models of the ther¬ 
mal Casimir force can be as large as one thousandt [40] . 
Preliminary results of an ongoing experiment at lUPUI 
based on the scheme of m, already prove conclusively 
that in the case of Ni the observed signal is three or¬ 
ders of magnitude smaller than the theoretical predic¬ 
tion based on the Drude model. In this paper we show 
how isoelectronic setups can be designed to further in¬ 
vestigate whether magnetic properties affect at all the 
Casimir force. 

The plan of the paper is as follows. In Sec. II we 
review the Drude and the plasma prescriptions for com¬ 
puting the thermal Casimir force. In Sec. Ill we present 
the general structure of isoelectronic setups for Casimir 
experiments, while in Sec. III-A and III-B we describe 
in detail setups specifically designed to investigate the 
thermal Casimir effect for non-magnetic and for mag¬ 
netic metals, respectively. Finally Sec. IV presents our 
conclusions. 


II. DRUDE AND PLASMA PRESCRIPTIONS 


According to Lifshitz theory [5], the Casimir free- 
energy T (per unit area) between two dieletric plane- 
parallel slabs Sj, j = 1,2 at distance a in vacuum is 
given by the formula: 
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where ks is Boltzmann constant, = 2TTlkBT/h are 
the (imaginary) Matsubara frequencies, k± is the modu¬ 
lus of the in-plane wave-vector, qi = \/, and 
k±) is the familiar Fresnel reflection coefficient 
of slab j for polarization a: 
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where = \/ej(i ^)^i /c^ + k']_, and Cj is the (dynamic) 
electric permittivity of slab j. According to Lifshitz 
formula, to compute the Casimir energy one needs to 
know the permittivty e(i^) of the involved materials along 
the imaginary frequency axis. This quantity cannot of 
course be measured directly, but it can be computed us¬ 
ing dispersion relations on the basis of optical data refer¬ 
ring to real frequencies w. For the case of insulators or 
ohmic conductors, the standard formula is provided by 
the Kramers-Kronig dispersion relation which expresses 
e(i^) in terms of the imaginary part e"(w) of the permit¬ 
tivity: 


e(ie) = l+ ^ 

TT 
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We see from Eq. (HI that in order to evaluate e(i^) at any 
imaginary frequency ^ it is in principle necessary to know 
e"(a;) at all frequencies w. Unfortunately such a com¬ 
plete knowledge of e"(w) is never possible, because optical 
data are always restricted to some finite frequency range 
Wmin < to < Wmax) Starting from a non-zero minimum fre¬ 
quency Wmin >0. In practice, there is no real difficulty 
on the high frequency side, because fall-off properties of 
e"(a;) ensure that for the relevant ^’s real frequencies w 
larger than a few tens of eV/fi give already a negligi¬ 
ble contribution to the integral on the r.h.s. of Eq. ([^. 
On the low-frequency side, however, one faces a problem 
when metals are considered. Since the imaginary part of 
the permittivity of ohmic conductors diverges like 1 /w 
for small frequencies, the integral in Eq. § receives a 
large contribution from low frequencies for which no op¬ 
tical data are available. Since truncation of the integral 
to the frequency Wmin would result in a large error, one 
is forced to extrapolate the dielectric function to 

frequencies uj < Wmin, where optical data are not avail¬ 
able, to evaluate the integral for ui < Wmin- As a rule, 
the extrapolation is done using the simple Drude model 
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where ujp is the plasma frequency, and 7 is the relaxation 
frequency. According to the Drude prescription, e(i^) is 


then estimated by the formula 

we"(w) 
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where £((,^.( 0 ;) is the imaginary part of the Drude permit¬ 
tivity Eq. (IM. The error introduced by this extrapola¬ 
tion in the estimate of e(i^), and thereof of the Casimir 
force, has been a subject of intense discussion [46l l47]. 
The problem can be partly relieved by using weighted 
dispersion relations [IHl SS], which sensibly reduce the 
contribution of the extrapolation. 

Surprisingly the estimate of e(i^) obtained by the above 
procedure, based on a Drude extrapolation of the opti¬ 
cal data, when plugged into Lifshitz formula results in a 
prediction of the Casimir force that appears to be incon¬ 
sistent with several recent experiments [IIHIIIIH]. It has 
been claimed that a prediction of the Casimir force which 
is consistent with the data can be obtained if e(i^) is com¬ 
puted by a different procedure, in which relaxation prop¬ 
erties of conduction electrons are neglected altogether. 
According to this so-called plasma prescription e(i^) is 
computed by the formula: 
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The main difference between the two prescriptions is in 
the power of divergence of £(iC) for vanishing While 
with the Drude prescription £(i^) has a simple pole, char¬ 
acteristic of ohmic conductors: 

^2 

^(iC) = —7 + 0(^°) (Drude prescription), (9) 
7? 


the plasma prescription leads to a double pole, similarly 
to superconductors m- 


^(iC) = ^ + (plasma prescription) . ( 10 ) 
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When the two prescriptions are used to compute the 
Casimir free energy between two metallic plates for room 
temperature, one gets almost coinciding values for the 
contributions of the non-vanishing Matsubara terms (i.e. 
the terms with Z ^ 0 in Eq. ([^). The I = 0 term with 
TM polarization is of course the same, because for both 
prescriptions the zero-frequency reflection coefficient for 
TM modes is one. The major difference is seen in the 
I = 0 TE mode, since this mode vanishes according to 
the Drude prescription while it does not according to the 
plasma prescription, as a result of the different singular 
behaviors of £(i^) for vanishing To give the reader 
a feeling of the magnitude of the difference among the 
Drude and plasma predictions, in Fig. Q we plot the 
corresponding reduction factors ij = J^/Eid for two gold 
plates, where Aid = —n^hc/{720a^) is the T = 0 energy 
(per unit area) between two perfectly conducting paral¬ 
lel plates. The solid line in Fig. Q is computed using 
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FIG. 1: Reduction factor rj — T /for two gold plates at 
room temperature, versus plates separation (in micron). The 
solid line is computed using the Drude prescription, while the 
dashed line is for the plasma prescription. The dot-dashed 
line corresponds to the plasma prescription, after subtraction 
of the contribution of the TE I = 0 mode. 

the Drude prescription, while the dashed line is for the 
plasma prescription. The dot-dashed line corresponds to 
the plasma prescription, after subtraction of the contri¬ 
bution of the TE I = 0 mode. The solid line is barely 
distinguishable from the dot-dashed line, which shows as 
anticipated earlier that the difference between the Drude 
and the plasma prescriptions arises solely from the ther¬ 
mal TE I = 0 mode. The plot shows also that the largest 
difference between the two prescriptions is found for sep¬ 
arations larger than one micron, where the thermal force 
becomes dominant. Compared to the plasma model, the 
Drude prescription predicts an extra repulsive thermal 
force which, was has been shown to originate from the 
magnetic interaction among thermal Eocault currents ex¬ 
isting in metallic plates an non-zero temperature mm- 
The above picture of the Drude vs. plasma compar¬ 
ison, described here for the case of two non-magnetic 
dielectric slabs, remains qualitatively true in the more 
general situations considered later in this work, when 
we move on to layered slabs and/or magnetic materi¬ 
als: in all cases the different magnitudes of the thermal 
force predicted by the Drude and the plasma prescrip¬ 
tions originate from the TE I = 0 mode, because the 
two prescriptions entail markedtly different values for the 
reflection coefficient of the TE mode for vanishing fre¬ 
quency. 


III. ISOELECTRONIC CASIMIR SETUPS 

As we discussed earlier, patch potentials constitute a 
major obtacle towards the observation of the thermal 
Casimir force. In recent searches of non-newtonian grav¬ 
ity [42l |43] the lUPUI group developed an elegant exper¬ 
imental technique based on isoelectronic force-difference 
measurements, which is by design immune from the prob¬ 
lem of potential patches. The structure of the isoelec¬ 


tronic setups we consider in this work is schematically 
illustrated in Fig. It consists of a sphere of radius 
R covered by a coating of material A and a planar slab 
divided in two regions, made of two different materials 
B and C, respectively. The key feature of the isoelec¬ 
tronic apparatus is the conductive over-layer of thickness 
w made of material D, covering both the B and the C 
regions. For any fixed sphere-plate separation a, we con¬ 
sider measuring the difference 

AF{a) = FB{a) - Fc{a) (11) 

among the values Fb and Fc of the (normal) Casimir 
force that obtain when the tip of the sphere is respec¬ 
tively above a point q deep in the B region, and a point p 
deep in the C region. The great advantage of this differ¬ 
ential mesurement over an absolute force measurement is 
that the detrimental (mean) electrostatic force caused by 
patches on the exposed surfaces of the plates is automat¬ 
ically subtracted out from AF, provided of course that 
the surface of the over-layer has identical patch struc¬ 
ture above the two regions of the plate. The limitations 
on the sensitivity of isolelectronic experiments resulting 
from random spatial variations of patch forces have been 
recently studied in detail [H], confirming the high sup¬ 
pression of patch forces in isoelectronic setups. An im¬ 
portant difference between our setups and those used by 
the lUPUI group should be stressed, however. In order to 
observe small differences between the gravitational inter¬ 
action of the sphere with the B and C sectors of the plate, 
the conductive overlayer of |121 05] was designed to be 
opaque, to screen out altogether the otherwise dominant 
electrostatic and Casimir forces. The thickness {w = 150 
nm) of the Au overlayers used in those experiments was 
therefore chosen to be larger than the plasma length of 
Au (Ap=130 nm). In our setup, instead, we do want 
to observe the differential thermal Casimir interaction of 
the sphere with the B and C regions of the plate. The 
conducting over-layer should then be designed such as 
to screen out the unwanted electrostatic component of 
the force, but at the same time it should impede as lit¬ 
tle as possible the passage of thermal photons that are 
responsible of the thermal Casimir interaction with the 
underlying materials. To do that, our overlayer has to 
be semi-transparent to infra-red radiation, and therefore 
its thickness has to be smaller than the skin depth of 
thermal photons for material D. 

We assume, as it is usually the case in Casimir ex¬ 
periments, that the sphere radius R is much larger than 
the separation a. Under this condition, it is possible to 
estimate the Casimir force by the Proximity Force Ap¬ 
proximation (PEA). The PEA has been widely used to 
interpret most Casimir experiments [6] (see |53j for more 
applications of the Proximity Approximation). Recently, 
it has been shown that the PEA represents the leading 
term in a gradient expansion of the Casimir force, in pow¬ 
ers of the slopes of the bounding surfaces [54H56] . Very 
recently the gradient expansion has been used to compute 
curvature corrections to the Casimir-Polder potential be- 
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FIG. 2: The isoelectronic setup consists of a sphere coated 
with material A above a planar slab divided in two regions, 
respectively made of materials B and C. Both the B and C re¬ 
gions are covered with a semi-transparent plane-parallel con¬ 
ducting over-layer of uniform thickness w made of material 
D. 


tween a small polarizable body and a curved surface m- 
The gradient expansion shows that the PFA is asymptot¬ 
ically exact in the zero-curvature limit, and with its help 
it is now possible to estimate the error introduced by 
the PFA. Besides R ^ a, two further conditions are re¬ 
quired to ensure the validity of the PFA in our setup. 
To be definite, we let (x, y, z) be cartesian coordinates 
such that {x, y) span the exposed surface of the conduc¬ 
tive overlayer, placed at z = 0, while the sphere tip is 
at z = a. We imagine that the a; < 0 region of the slab 
is made of material B, while its a; > 0 region is made of 
material C. To ensure that we can neglect the effect of 
the sharp boundary between the B and the C regions, 
we assume that the horizontal distances of the points p 
and q from the B-C boundary are both much larger than 
the typical interaction radius p = V aR of the circular 
region of the plate that contributes significantly to the 
Casimir interaction: s ^ p. The force Fb (Fc) is then 
undistinguishable from the force Fb {Fq) between the 
sphere and a plane-parallel two-layers slab consisting of 
a layer of thickness w of material D on top of an infi¬ 
nite dielectric planar slab made of material B (C). Under 
these conditions, we have for AF: 

AF ~ Ub - Fc ~ 


= 2ttR{Fb - Fc) , (12) 

where and denote the PFA expressions 

for the Casimir force between the sphere and a two-layer 
slab consisting of the conductive layer of material D on 
top of either a B or a C substrate, respectively, while Fb 
and Fc denote the unit-area Casimir free energies for the 


corresponding plane-parallel systems in which the sphere 
is replaced by a planar slab of material A. In the last 
passage of Eq. (12) we used the well-known PFA for¬ 
mula for the Casimir free energy of a sphere-plate system 
— 27rFF(a). For later applications, we con¬ 
sider the possibility that some among the materials A, B 
and C are magneto-dielectric. In its original version Lif- 
shitz theory was formulated for planar dielectrics {p = 1) 
fully described by the respective frequency dependent 
(complex) dynamical permittivity e(a;). The theory was 
later generalized to deal with layered magneto-dielectric 
plates in [^. It turns out that the Casimir free- 
energy F (per unit area) between two magnetodieletric 
possibly layered plane-parallel slabs Sj, j = 1,2 at dis¬ 
tance a in vacuum is still represented by Eq. pro¬ 
vided that the reflection coefficients are now understood 
to be those for the possibly layered slabs. The Casimir 
free energy Fb can be obtained from Eq. (§ by substitut¬ 
ing Ra'^ by the Fresnel reflection coefficient (given 

in Eqs.(14) and (151 below, with a = 0, b =A), and Ra^ 
by the reflection coefficient of a two-layer planar 

slab consisting of a layer of thickness w of material D on 
a dielectric B substrate. The latter reflection coefficient 
has the expression: 


Fi°°^)(iCi,fc_L;u;) 


,PB) 
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where = \/e.a(p^i)+ k^ , ta and pa de¬ 
note the electric and magnetic permittivities of medium 
a, and we define cq = pp = 1. The Casimir free energy 
Fc for a planar A-D-C system can be obtained by re¬ 
placing material B by material C in the above formulae. 


A. Non-magnetic metals 

The first setup we consider is designed to discriminate 
among the Drude and the plasma prescriptions for non¬ 
magnetic conductors. For this purpose we found that a 
convenient choice of materials is the following: we take 
Au for the coating of the sphere and for the C region 
of the plate, and high-resistivity (dielectric) Si for the 
region B. The choice of the material for the overlayer is 
critical. As we explained earlier, the over layer needs to be 
semitransparent to thermal photons. This requirement 
led us to consider B-doped low-resistivity Si, because its 
large plasma length (Ap = 2.7^m) ensures a good degree 
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FIG. 3: Force difference AF (in fN) versus separation a 
(in microns) for a An sphere of radius R = 150 pm in front 
of a Si-Au plate, covered with 100 nm of low-resistivity Si. 
The red and grey lines correspond to the plasma and Drude 
prescriptions, respectively. 


of transparency also for relatively large thicknesses w (we 
take w = 100 nm). B-doped Si has been used successfully 
in Casimir experiments (see the second of Refs. |15)). 
Below, we shall use the symbol Sic to denote conductive 
silicon, while the symbol Si shall denote dielectric silicon. 

The values €a{i^i), (a=Au, Si) of the permittivities of 
Au and dielectric Si were computed by the procedure 
expalined in Sec. II, on the basis of the tabulated optical 
data quoted in [5D]. When using the Drude prescription, 
the permittivity eAu(iC) of Au was computed according 
to Eq. Q. The optical data for Au were extrapolated 
towards low frequencies by the Drude model, with Drude 
parameters wau = 8.9eV/h, 7 au = 0.035 eV/7i. For the 
permittivity of conductive Si we used the formula 


esic(iC) = esi(iC) + t > 

?(? + 7sicJ 


(16) 


with plasma frequency wsp = 0.46 eV /h and relaxation 
frequency 7 si„ = O.leV/S (see Ref. [5], pag. 588). When 
we considered the plasma prescription, we made use of 
Eq. ([^ to compute the imaginary-frequency permittivity 
of Au. In the case of Sic the plasma prescription was 
implemented simply by setting to zero 731 ,, in Eq. (16) . 

In Fig. 1^ we plot AF (in fN) versus separation a 
(in /rm) for a sphere of radius i? = 150 p,m and for a 
thickness w = 100 nm of the Sic over-layer. The red 
and the grey curves correspond the plasma and Drude 
prescriptions, respectively. The plasma and the Drude 
models predict widely different magnitudes for AF. For 
example, for a = 3 ^m AFpiasma is fourteen times larger 
than AEDrude, while for a = 4 /rm they differ by a fac¬ 
tor around fifty. This shows that our setup produces 
a strong amplification of the Drude-plasma discrepancy, 
which makes one confident that an unambiguous discrim¬ 
ination between the two models should be possible with 
this apparatus. The recent isoelectronic differential ex¬ 
periment at lUPUI [33], using a Au coated sphere of ra¬ 
dius R ~ 150/im, achieved a sensitivity better than 0.3 


FIG. 4: Force difference AF (pN) versus separation a (in 
nm) for a Ni-coated sphere and a Au-Ni plate, covered by 80 
nm of Au. The curve was computed using the Drude model 
with pNi(O) = 110. 


AF(10“‘^ N) 



a (nm) 


FIG. 5: Force difference AF (in fN) versus separation a 
(in nm) for a Ni-coated sphere and a Au-Ni plate, covered 
by 80 nm of Au. The black line corresponds to the Drude 
model with ^Ni(O) = 1, while the solid and dashed red lines 
correspond to the plasma model, with /rNi(O) = 110 and 
MNi(O) = 1, respectively. 


fN in force differences, independent of the separation a in 
the range from 200 nm to 1 /tm. With this level of sen¬ 
sitivity, it would be possible to accurately measure AF 
up to separations of several ^m. 


B. Magnetic metals 

In this Section we consider setups involving a ferro¬ 
magnetic material. Magnetic materials bring up an in¬ 
teresting problem in connection with the Casimir effect: 
does the large magnetic permeability of ferromagnetic 
substances affect the Casimir force? To address this ques¬ 
tion, we recall that the dynamic magnetic permittivity of 
magnetic materials can be modeled by the Debye formula 


m 


/r(a;) — 1 -|- 


m -1 

1 - iuj/uJrr 


(17) 
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AF(10“‘^ N) 



FIG. 6: Force difference AF (in IN) versus separation a (in 
nm) for a Ni-coated sphere and a Pt-Ni plate, covered with 20 
nm of Pt. The black line corresponds to the Drude model with 
^i,ji(0) = 1, while the solid and dashed red lines correspond 
to the plasma model, with /rNi(O) = 110 and /iNi(O) = 1, 
respectively. 
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FIG. 7: Force difference AF (in fN) versus separation a (in 
microns) for a Ni-coated sphere and a Pt-Ni plate, covered 
with 100 nm of conductive Si. The black line corresponds to 
the Drude model with ^Ni(O) = 1, while the solid and dashed 
red lines correspond to the plasma model, with pNi(O) = 110 
and pNi(O) = 1, respectively. 

AF(10“‘^ N) 



FIG. 8: Force difference AT (in fN) versus separation a (in 
microns) for a Ni-coated sphere and a Si-Ni plate, covered 
with 100 nm of conductive Si. The black line corresponds to 
the Drude model with pNi(O) = 1, while the solid and dashed 
lines correspond to the plasma model, with /rNi(O) = 110 and 
MNi(O) = 1, respectively. 


where fj,(0) is the static magnetic permeability and ojm 
is a characteristic frequency, typically smaller than 10® 
Hz. Since is several orders of magnitude smaller than 
the frequency of the first Matsubara mode ^ 10^^ Hz 
at room temperature, when computing the Casimir free 
energy one can set /i = 1 in alH > 0 Matsubara terms 
in Eq. (H). This means that magnetic properties can 
affect the Casimir force only via the I = 0 TE term in 
Eq. (§. The latter term is proportional to the temper¬ 
ature, and therefore it is clear that this term must de¬ 
scribe the effect of thermal fluctuations of the magnetic 
dipoles inside the ferromagnetic plates. However, it has 
been known for a long time that magnetic domains ex¬ 
isting within (non magnetized) magnetic materials show 
thermal fluctuations only in small particles or at tem¬ 
peratures close to the Curie point Tc, while under nor¬ 
mal circumstances the equilibrium magnetic microstruc¬ 
ture is athermal [62] . This makes one wonder whether 
magnetic properties affect the Casimir force at all, far 
from the Curie point. Below we shall consider specifi¬ 
cally the case of Ni. This ferromagnetic material has a 
static permeability /iNi(O) = 110, and a Curie tempera¬ 
ture Tc = 620 K. In order to address the question of the 
influence of the static permeability of Ni on the Casimir 
force, we shall consider within both the Drude and the 
plasma prescriptions two further alternatives to estimate 
the I = 0 TE mode, i.e. to take /rNi(O) = 110 or instead 
to take /rNi(O) = 1, the two alternatives corresponding re¬ 
spectively to full inclusion or full neglect of the magnetic 
properties in the Casimir force. 

Recent experiments by the Riverside group measured 
the gradient of the Casimir between a Au sphere and a Ni 
plate and between a Ni coated sphere and a Ni plate 
[IS] , in the separation range from 200 nm to 550 nm. The 
data of these experiments were interpreted as providing 
evidence for the influence of the magnetic properties of 
Ni on the Casimir force, and to be in agreement with the 
plasma prescription but not with the Drude model. As 
it is usual in Casimir experiments operating in the sub¬ 
micron range, the relative difference between the plasma 
and the Drude models was found to be rather small, less 
than about ten percent. Therefore, the theoretical in¬ 
terpretation of the data required a delicate analysis of 
several possible sources of error. 

We show in this Section that isoelectronic setups 
strongly enhance the discrepancy between the Drude and 
the plasma models in magnetic systems, both for submi¬ 
cron separations and for micron separations of the plates. 
We shall consider several such setups, with the purpose 
of obtaining a complete discrimination between the four 
different possible theoretical models that were considered 
in [ndS], i.e. the Drude and the plasma models with 
account or neglect of the magnetic properties of Ni. 

The first setup is designed to discriminate the Drude 
model with p,Ni(0) = 110 from the three other models: it 
uses Ni both for the sphere coating and the region C of 
the plate, while the overlayer and region B are both made 
of Au. The permittivity eAu(iC) of Au was computed by 















the same procedure described in the previous Section. 
The permittivity eNi(iC) of Ni was computed by a similar 
procedure, using the tabulated optical data of m ex¬ 
trapolated to low frequencies via a Drude model with pa¬ 
rameters WNi = 4.89 eV/h and 7 Ni = 0.0436 eV/H [TS]. In 
Figs. I^andj^we plot AF versus separation a (in nm), for 
a sphere of radius i? = 150 /im and for a thickness ic = 80 
nm of the Au overlayer. The curve in Fig. was com¬ 
puted using the Drude prescription with /rNi(O) = 110, 
while the three curves in Fig. were computed using 
the Drude prescription with /rNi(O) = 1 (black line), and 
the plasma prescription with /rNi(O) = 110 (red line) and 
MNi(O) = 1 (dashed red line). Remarkably, we see that 
the Drude model with ^Ni(O) = 110 predicts a signal 
whose magnitude is three orders of magnitude larger than 
that predicted by the three other models. Therefore, by 
this setup is should be easy to discriminate the Drude 
model with with /rNi(O) = 110 from the other three mod¬ 
els. An experiment based on this scheme is presently on¬ 
going at lUPUI. The setup is identical to that described 
in [43], apart from the replacement of the gold coated 
sphere by a Ni coated one, and the replacement in the 
rotating disk of the Si sectors by Ni sectors. Preliminary 
data presented by R.S. Decca at a conference in Cancun 
(Mexico) show that for a sample using a 84 nm Au over¬ 
layer AF is of the order of -1 fN (with an error of 0.3. fN) 
in the range of separations from 200 nm to 400 nm. The 
observed AF has the opposite sign and is off by three or¬ 
ders of magnitudes with respect to the prediction of the 
Drude model with /rNi(O) = 110, while it is in qualitative 
agreement with the remaining three models. This rules 
out conclusively the Drude model with full account of the 
Ni permeability. It remains to see if the Ni permeability 
has any impact at all on the Casimir force. To achieve a 
good discrimination between the remaining three models 
that we consider, it is necessary to change some of the 
materials used in the first setup. We consider first a setup 
in which Au is replaced by Pt, both for the overlayer and 
for region B. The permittivity ept(iC) of Pt was computed 
by the same procedures described earlier, using the tabu¬ 
lated optical data of m extrapolated to low frequencies 
via a Drude model with parameters wpt = 5.1eV/7i and 
7Pt = 0.07eV/7i [S3]. In Fig. [^we show a plot of AF 
(in FN) versus separation (in nm) for a sphere of radius 
R = 150 ^m, and for a thickness w = 20 nm of the Pt 
overlayer. We see that in the range from 200 to 300 nm 
the plasma model with ^Ni(O) = 110 predicts a positive 
AF, while both the Drude and the plasma models with 
AiNi(O) = 1 predict a negative AF. Such a sign difference 
should be easy to detect. Observation of a positive AF 
would provide a clear signature of the influence of the 
magnetic permeability on the Casimir force. Another 
setup that could allow for a clear discrimination of the 
plasma model with /rNi(O) = 110 versus the two other 
models with /rNi(O) = 1 uses a 100 nm overlayer made 
of conductive Si, and again a region B made of Pt. In 
Fig. 0 we show a plot of AF (in FN) versus separation 
(in micron) for a sphere of radius R = 150 fim. We see 


from the plot that with this setup also it should be easily 
possible to discriminate the /rNi(O) = 110 plasma model 
(solid red line) from the two /rNi(O) = 1 models (black 
and dashed red line). A good discrimination among the 
latter two models is finally achieved by the third setup, 
which uses a Ni coated sphere in front of a plate made 
of dielectric Si (B region) and Ni (C region), covered by 
a conductive Si overlayer. In Fig. |^we show a plot of 
AF (in FN) versus separation (in micron) for a sphere 
of radius R = 150 p,m and for a thickness w = 100 nm 
of the Sic overlayer. We see that with this setup there is 
indeed a large difference between the signals predicted by 
AiNi(O) = 1 Drude model (grey curve) and the ^Ni(O) = 1 
plasma model (dashed red-line). We conclude that the 
combined use of the three Ni setups should allow to de¬ 
termine whether or not the magnetic permeability of fer¬ 
romagnetic substances has any infuence on the Casimir 
effect, and which among the plasma and Drude models 
better describes it. 


IV. CONCLUSIONS 

The problem of thermal corrections in the Casimir ef¬ 
fect still constitutes a major unresolved puzzle in this 
field. Different theoretical prescriptions have been pro¬ 
posed to compute the thermal Casimir force, which pro¬ 
vide different predictions for the thermal force. Several 
experiments have been performed in recent years to clar¬ 
ify the matter, but the situation remains unclear. Some 
experiments [nHHUH] probing the Casimir force for sep¬ 
arations between 200 and 700 nm have been interpreted 
as supporting the plasma prescription, and to rule out 
the Drude prescription, while a single experiment |38j 
probing the wide range from 700 nm to 7.3 micron was 
interpreted to be consistent with the Drude model, and 
inconsistent with the plasma model. A major problem 
in the theoretical interpretation of Casimir experiments 
is represented by stray electrostatic forces caused by po¬ 
tential patches on the plates. This represents a severe 
problem, especially for micron separations of the plates 
where stray forces become large. 

An experimental approach which is in principle im¬ 
mune from the problem of patch potentials is based on 
isoelectronic force difference measurements. This ap¬ 
proach has been successfully used by the lUPUI group in 
recent searches of non-newtonian gravitational forces in 
the submicron range [42l|43|. The author showed recently 
snun] that the isoelectronic technique can be adapted 
to investigate the thermal Casimir effect. In this paper 
we have further developed this idea. We have shown 
that by suitably choosing the materials of the samples, 
it is possible to investigate in great detail the features of 
the thermal Casimir force, both for magnetic and non¬ 
magnetic metals. Apart from being immune from elec¬ 
trostatic problems, the devices described in this paper 
provide a strong amplification of the differences between 
the Drude and the plasma prescriptions, thus making 
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a discrimination among these models much simpler and 
hopefully conclusive. 

Preliminary results of an ongoing experiment at lUPUI 
with Ni samples [IS], based on the isoelectronic scheme 
described in this paper, already prove conclusively that 
the effect of the magnetic permeability on the Casimir 
force is three orders of magnitude smaller than what can 
be predicted on the basis of the Drude model. The avail¬ 
able data are qualitatively consistent with the plasma 
model. Another possibility to interpret the data is that 
magnetic propertites do not affect at all the Casimir ef¬ 
fect, in which case both the Drude and the plasma models 
are qualitatively in agreement with the data. In this pa¬ 
per we have described three different isolectronic setups 
that should allow to fully resolve this problem, and to es¬ 


tablish which among the Drude and the plasma models 
provides the correct description of the Casimir effect in 
magnetic systems. 
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